Fermionic bound states on a one-dimensional lattice 
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We study bound states of two fermions with opposite spins in an extended Hubbard chain. The 
particles interact when located both on a site or on adjacent sites. We find three different types 
of bound states. Type U is predominantly formed of basis states with both fermions on the same 
site, while two states of type V originate from both fermions occupying neighbouring sites. Type 
U, and one of the states from type V, are symmetric with respect to spin flips. The remaining one 
from type V is antisymmetric. V-states are characterized by a diverging localization length below 
some critical wave number. All bound states become compact for wave numbers at the edge of the 
Brilloin zone. 

PACS numbers: 03.65.Ge, 34.50.Cx, 37.10.Jk 
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Introduction. Advances in experimental techniques of 
manipulation of ultracold atoms in optical lattices make 
it feasible to explore the physics of few-body interactions. 
Systems with few quantum particles on lattices have new 
unexpected features as compared to the condensed mat- 
ter case of many-body interactions, where excitation en- 
ergies are typically small compared to the Fermi energy. 
In particular, a recent experiment explored the repulsive 
binding of bosonic atom pairs in an optical lattice [l| , as 
predicted theoretically decades earlier 0, [1] (see also [1| 
for a review). In this work, we study binding properties of 
fermionic pairs with total spin zero. We use the extended 
Hubbard model, which contains two interaction scales 
- the on site interaction U and the nearest neighbour 
intersite interaction V. The nonlocal interaction V is 
added in condensed matter physics to emulate remnants 
of the Coulomb interaction due to non-perfect screening 
of electronic charges. For fermionic ultracold atoms or 
molecules with magnetic or electric dipole-dipole inter- 
actions, it can be tuned with respect to the local interac- 
tion U by modifying the trap geometry of a condensate, 
additional external dc electric fields, combinations with 
fast rotating external fields, etc (for a review and relevant 
references see Q). 

The paper is organized as follows. We first describe 
the model and introduce the basis we use to write down 
the Hamiltonian matrix to be diagonalized. Then we de- 
rive the quantum states of the lattice containing one and 
two fermions, with opposite spins in the latter case. We 
study the obtained bound states for two fermions. We 
obtain analytical expressions for the energy spectrum of 
the bound states. Some of the bound states are charac- 
terized by a critical momentum below which they dissolve 
with the two-particle continuum. 

Model and basis choice. We consider a one-dimensional 
lattice with / sites and periodic boundary conditions de- 
scribed by the extended Hubbard model with the follow- 
ing Hamiltonian: 



where 



a 



Hi 



(2) 
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Hq describes the nearest-neighbor hopping of fermions 
along the lattice. Here the symbols a —1,1 stand for a 
fermion with spin up or down. Hu describes the onsite 
interaction between the particles, and Hy the intersite 
interaction of fermions located at adjacent sites, dj'^ and 
are the fermionic creation and annihilation operators 



a 



satisfying the corresponding anticommutation relations: 
^+ 1 - -^--^ ■ %+^,a^<,,] = [dj,a.di..'\ = 0. Note 



H = Ho + Ht 



Hi 



(1) 



that throughout this work we consider U and V posi- 
tive, which leads to bound states located below the two- 
particle continuum. A change of the sign of U, V will 
simply swap the energies. 

The Hamiltonian ([T]) commutes with the number op- 
erator N = 'Y^- hj whose eigenvalues are n = + n^, 
i.e. the total number of fermions in the lattice. We con- 
sider n — 2, with rij = 1 and = 1. Therefore we 
construct a basis starting with the eigenstates of TV. We 
use a number state basis = |ni; n2 • ■ ■ n/) Qj where 
TT-i — Ui^l + rii^i represents the number of fermions at the 
i-th site of the lattice. \^n) is an eigenstate of the number 
operator N with eigenvalue n = X]J=i ^j- 

To observe the fermionic character of the considered 
states, any two-particle number state is generated from 
the vacuum \0) by first creating a particle with spin 
down, and then a particle with spin up: e.g. aJ|d^JO) 
creates a particle with spin down on site 1 and one with 



spin up on site 2, while 0,2 ^\0) creates both particles 
with spin down and up on site 2. 

Due to periodic boundary conditions the Hamilto- 
nian ([T]) commutes also with the translation operator T, 
which shifts all lattice indices by one. It has eigenval- 
ues r = exp{ik), with Bloch wave number k = and 
1/ = 0,1,2,...,/- 1. 

Single particle spectrum. For the case of having only 
one fermion (either spin up or spin down) in the lattice 
(n = 1), a number state has the form \j) — a^^|0). The 

interaction terms Hu and Hy do not contribute. For a 
given wave number fc, the eigenstate to ([!]) is therefore 
given by: 

^ s— 1 

The corresponding eigenenergy 

Efc — —2 cos(fc). (6) 

Two fermions with opposite spins. For two particles, 
the number state method involves Ng = basis states, 
which is the number of ways one can distribute two 
fermions with opposite spins over the / sites including 
possible double occupamcy of a site. Below we consider 
only cases of odd / for simplicity. Extension to even 
values of / is straightforward. 

We define basis states to a given value of the wave 
number k: 

/ 

l*i) = ;^E(7)'"'«tT«ulo>, (7) 

\''^^+)=4ji:{^y"-tA\o), (8) 

^ s— 1 

/ 

\'^^^-)-^i:{^y"-tAi\o), (9) 

l<J<(/ + l)/2. (10) 

Note, that the sign ± discriminates between states, where 
the distance from the spin up particle to the spin down 
particle is smaller (larger) than vice versa. Note that 
distances are measured by scanning the (periodic) chain 
in the direction of increasing lattice site number, passing 
a given particle (say with spin up) and then counting the 
distance to the spin down particle. In the limit of an 
infinite system the two signs discriminate between states 
where the spin up particle is to the left or right of the 
spin down one. 

Therefore, a complete wavefunction is given by 

f±l f+i 

2 2 

1*2) = cii$i) + J2 + E ■ (11) 

Any vector in our given Hilbert space is then spanned by 
the numbers |ci, C2+, C2-, C3+, ...). 
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FIG. 1: Energy spectrum of the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. 
Here f/ = l, V = 1, / = 101. Symbols are from analytical 
derivation, lines are the result of numerical diagonalization. 
The arrows indicate the location of the critical wave numbers 
(see text). 

Next we calculate the matrix elements of the Hamilto- 
nian ^ in the framework of the basis pT|) . We arrive at 
a / X / matrix with elements H(i,j) (i, j — I, . . . , d): 

/U q q ... \ 

q* V q ... 

q* V q ... 

9* g 



... q* q 
... q* Op 
\ ... q* p Oj 

Here g = 1 + r and p = r^^^+i^/^ + r~(/-i)/2. 

In Fig. [T]we show the energy spectrum of the Hamilto- 
nian matrix (|12p obtained by numerical diagonalization 
for the interaction parameters U — 2 and I^ = 2 and 
/ = 101. At J7 = and V ~ 0, the spectrum is given by 
the two fermion continuum, whose eigenstates are charac- 
terized by the two fermions independently moving along 
the lattice. In this case the eigenenergies are the sum of 
the two single-particle energies: 

<,fc2 = -2[cos(/ci)-f cos(/c2)], (13) 

with fci,2 = 7rj/i,2/(/ + 1) , i^i,2 - 1, . . . , /. The Bloch 
wave number k = ki+k2 mod 2n. Therefore, if k — ±7r, 
the continuum degenerates into points. The continuum is 
bounded by the hull curves h±{k) — ±4 cos |. The same 
two-particle continuum is still observed in Figure[T] for 
nonzero interaction. However, in addition to the contin- 
uum, we observe one, two or three bound states dropping 
out of the continuum, which depends on the wave num- 
ber. For any nonzero U and V, all three bound states 
drop out of the continuum at fc = ±tt. One of them stays 
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A Antisymmetric bound states 
□ Symmetric bound states 
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FIG. 2: Energy spectrum of the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. 
Here U = 4 , V = 3, f — 101. The symbols are from 
analytical derivation, lines are the results of numerical diag- 
onalization. The arrow indicates the location of the critical 
wave number (see text). 



bounded for all values of k. The two other ones merge 
with the continuum at some critical value of upon 
approaching fc = as observed in FiglT] Note that for 
k = ±7r and U = V, all three bound states are degener- 
ate. 

Upon increasing U and V, we observe that a second 
bound state band separates from the continuum for all k 
(Figl21). At the same time, when U ^ V, the degeneracy 
at fc = ±7r is reduced to two. 

Finally, for even larger values of U and V, all three 
bound state bands completely separate from the contin- 
uum (Figini). 

Symmetric and antisymmetric state representation. In 
order to obtain analytical estimates on the properties 
of the observed bound states, we use the fact that the 
Hamiltonian for a two fermion state is invariant under 
flipping the spins of both particles. We define symmetric 
basis states 



1 



and antisymmetric states 



(l<i>.,+> + l'5„-)) 



(14) 



(15) 



Note that |$i) is a symmetric state as well. Then the ma- 
trix (fT2|) can be decomposed into irreducible symmetric 
and antisymmetric ones: 



( U qV2 

q* 



\ 



q 
q* pj 



(16) 



FIG. 3: Energy spectrum for the two fermion states. The 
eigenvalues are plotted as a function of the wave number k. 
Here U = 8 , ^ = 8, / = 101. The symbols are from 
analytical derivation, lines are the results of numerical diago- 
nalization. 



q' 



V 





q* 



q 

-pJ 



(17) 



The rank of the matrix iJ* is (/ + l)/2, while the rank 
of i?'^ is (/- l)/2. 

Antisymmetric bound states. The antisymmetric 
states exclude double occupation. Therefore the spec- 
trum is identical with the one of two spinless fermions 
Following the derivations in Q we find that the an- 
tisymmetric bound state, if it exists, has an energy 



E^ik) = -{V + 



V 



(18) 



This result is valid as long as the the bound state en- 
ergy stays outside of the continuum. The critical value 
of k at which validity is lost, is obtained by requesting 
\E^{k)\ = \h±{k)\. It follows V = 2cos(|). Therefore 
the antisymmetric bound state merges with the contin- 
uum at a critical wave number 



V 

fc° = 2 arccos( — ) , 



(19) 



For F = 1 it 



setting a critical length scale A° = 
follows fc^/TT « 0.667 (see FigH]). 

The equation (fT5|) is in excellent agreement with the 
numerical data in Figs. I1I2I3I (cf. open triangles). We 
also note, that the antisymmetric bound state is located 
between the two symmetric bound states, which we dis- 
cuss next. 

Symmetric bound states. A bound state can be 
searched for by assuming an unnormalized eigenvector 



4 



to (|16p of the form |c, 1, /i, /i^, /i^, ...) with = p < 1. 
We obtain 

Ec=Uc+ V2q* , 
E = V2qc + V + q*fi, (20) 

E=- + q*ii . 
/i 

It foUows that /I — pcxp(ik/2) and 

E^{k) = -2{p+-) cos k/2 . (21) 
P 

The parameter p satisfies a cubic equation 

ap^ + bp'^ + cp + d = (22) 

with the real coefficients a, b. c and d given by a — 
2Fcos(|), b = 4cos2(|) - J7V^, c = 2{U + V)cos{^), d = 
— 4cos^(|). An analytic solution to can be obtained, 
but is cumbersome to be presented here. We plot the 
results in Figs. I1I2I3I (cf. open circles and squares). We 
obtain excellent agreement. 

At the Brilloin zone edge k = ±7r the cubic equation 
([22|) is reduced to a quadratic one, and can be solved to 
obtain finally p — > and 

E'2^{k^±i:) = -U , Ef{k^±'K) = -V . (23) 

In particular we find for k — ±7r that E^ — i?2. In 
addition, if C/ — V, all three bound states degenerate at 
the zone edge. 

If y = 0, the cubic equation (|22|) is reduced to a 
quadratic one in the whole range of k and yields Q 

£;|^ (fc) - - Vt^2 + i6cos2(fc/2) . (24) 

Next we determine the critical value of k for which the 
bound state with energy E^ is joining the continuum. 



Since at this point p = 1, we solve (|22|) with respect to 
kc and find 

setting another critical length scale = |f . E.g. for 

U = V = 1 kl^ /tt k, 0.89, in excellent agreement with 
FigUl For [/ = 4 and F = 3 we find fcf /tt w 0.59 
confirming numerical results in Figl2l 

Conclusions. Two fermionic particles with opposite 
spin allow for three different types of bound states on a 
one-dimensional lattice with onsite U and nearest neigh- 
bour V interaction. Two of them are symmetric with 
respect to spin flips, and one is antisymmetric. The anti- 
symmetric bound state is characterized by a critical wave 
number separates wave numbers with bound states from 
wave numbers without. It follows from that this hap- 
pens for F < 2. For larger values of V the whole wave 
number space becomes available for antisymmetric bound 
states, similar to one of the symmetric bound states for 
any nonzero U. The second symmetric bound state also 
observes a critical wave number. It follows from (|25p that 
this happens for U < AV/{V — 2), while the whole wave 
number space becomes available otherwise. It could be 
a challenging task to observe these different phases with 
one, two, or three bound states experimentally, by tuning 
[/, V , and k. 

For higher lattice dimensions more nearest neighbours 
have to be taken into account, similar to an increase 
of the interaction range. In these cases, we expect 
consequently more bound states to appear. 
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